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Abstract 

We report a new measurement of the branching fraction of neutral B meson decays 
S°^aj t (1260)vr =F with af (1260) -^p ^ and ^^tt+tt - using 275 x 10 6 BB pairs collected by the 
Belle detector at the KEKB asymmetric-energy e + e~ collider. We measure the branching fraction 
B(B° -► 0^(1260)71-^) = (48.6 ± A.l(stat) ± 3.9(syst)) x 10" 6 . Using a relativistic Breit-Wigner 
parameterization, we measure the mass and width of the ai(1260) to be m ai = 1197 ± 34 MeV/c 2 
and r ai = 305 ±43 MeV/c 2 , respectively. 



PACS numbers: 13.25.HW, 11.30.Er, 12.15.Hh, 14.40.Nd 
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The Kobayashi and Maskawa (KM) model explains the source of CP violation in terms 
of a single complex phase in the quark mixing matrix [l[. Measurements of the CP violating 
asymmetry parameter sin 20i by the Belle and BaBar || collaborations have established 
CP violation in the neutral B meson system. Measurements of other CP violating asymme- 
try parameters provide important tests of the KM model. Decay modes where the b — > uud 
contribution is dominant, such as B° — > 7r + 7r~, B° — > pre, B° — > pp and B° — > o^tt can 
be used to measure the Cabibbo-Kobayashi-Maskawa angle 02 (also denoted a) |4|. Un- 
like B° — > 7T + 7r _ decay, B° — > afir^ decay is not a CP eigenstate, and four flavor-charge 
configurations (B°(B°) — > af tt t ) must be considered to measure the angle 02. 

The CLEO collaboration has quoted an upper limit of 49 x 10~ 5 at the 90% C.L. for 
the branching fraction of B° — > af (1260)7r =F |5[, while the DELPHI collaboration has given 
an upper limit of 28 x 10~ 5 for the branching fraction of B° — > 4ir at the 90% C.L. 
Recently the BaBar collaboration has reported a measured branching fraction of B(B° —* 
af (1260)7r T ) = (40.2±3.9(stat)±3.9(syst)) x 10~ 6 [7]. In some previous studies discrepancies 
have been found in the ai(1260) properties. Among these studies, the analyses involving 
hadronic events and r decays Q are important. Therefore, it is important to measure 
the branching fraction of this decay mode. A new precise measurement of the ai(1260) 
properties is also useful. 

We present in this paper a measurement of the branching fraction of B° af (1260)ti t 
with 0^(1260) — > p°7r ± and p° — > 7r + 7r~ [10]. In this analysis, we assume that the main 
contributions come from the B° — > af(12Q0)-K T and B° —> 0^(1320)71: T decays, and neglect 
any interferences between decay modes consisting of four charged pions in the final state. 
We also assume that ai(1260) decays to pn only. 

The analysis is based on a 253 fb _1 data sample containing 275 million BB pairs. The 
data were collected with the Belle detector at the KEKB asymmetric-energy e + e~ (3.5 
on 8 GeV) collider [12] operating at the T(4S*) resonance (y/s = 10.58 GeV) with a peak 
luminosity that exceeds 1.5 x 10 34 cm _2 s _1 . The Belle detector is a large- solid- angle magnetic 
spectrometer that consists of a silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cerenkov counters (ACC), a barrel- like arrangement 
of time-of-flight scintillation counters (TOF) and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals located inside a super-conducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside the coil is instrumented to detect K\ 
mesons and to identify muons. Two inner detector configurations were used. A 2.0 cm 
radius beampipe and a 3-layer silicon vertex detector were used for the first sample of 152 
million BB pairs, while the remaining 123 million BB pairs were recorded using a 1.5 cm 
radius beam pipe, a 4-layer silicon detector and a small-cell inner drift chamber [l3j . 

For the present analysis several Monte-Carlo (MC) samples (signal and backgrounds) are 
generated with EvtGen 0). The Belle detector response is simulated with a GEANT3-based 
program [l5| . 

We reconstruct the B° — > af (1260)7r T candidate using four charged tracks originating 
from the beam interaction region and each having a momentum transverse to the beam, p t , 
greater than 100 MeV/c. To identify K and ir mesons, we form a K(n) likelihood L^iL^ 
by combining information from the CDC (dE/dx), the TOF and the ACC. Discrimination 
between pions and kaons is achieved through the likelihood ratio R(K/ii) = Lk/{L^ + Lk)- 
Charged tracks with R(K/tc) < 0.4 are regarded as pions. Positively identified electrons and 
muons are rejected. 

In the event selection we reconstruct each p° candidate by combining a tt + and a 7r~ 
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such that the invariant mass of the pion pair (171^+^-) satisfies 0.55 GeV/c 2 < m^^- < 
1.15 GeV/c 2 . We reconstruct the af (1260) candidates by combining another charged pion 
with the reconstructed p°. The af (1260) candidates are then selected on the invariant mass: 
0.80 GeV/c 2 < m ai < 1.775 GeV/c 2 . 

Candidate B mesons are identified using three kinematic variables: the beam-energy- 
constrained mass, Mb c = J -E 2 eam — Vbi the energy difference, AE = E B — E heSLm , and the 
helicity angle of the ai(1260) meson. Here, Eb and ps are the reconstructed energy and 
momentum of the B candidate in the center of mass (cm.) frame, and -Ebeam is the beam 
energy in the cm. frame. The angle between the p° momentum vector and the direction 
opposite to the B° in the af (1260) rest frame is defined as the helicity angle (#hci)- To 
extract the signal yield we select events in the region 5.21 GeV/c 2 < M bc < 5.29 GeV/c 2 , 
—0.10 GeV < AE < 0.12 GeV and — 1.0<cos#hei< 1-0 (loose signal region). The momentum 
of the bachelor pion (the pion that comes directly from the decay of B° — > af (1260)7r T ) in 
the B° rest frame is required to satisfy 2.2 GeV/c < Pbach < 2.7 GeV/c. We find 39% of the 
reconstructed events have more than one B candidate after all selection requirements. We 
select the best B meson candidate that has the minimum \ 2 value in the vertex fit to the 
four charged pion tracks. 

We find 14% of signal events in the MC simulation are incorrectly reconstructed by in- 
cluding at least one charged track from the other B meson. The remaining 86% of signal 
events are reconstructed from the correct set of final state pions. However there remains the 
possibility of incorrectly assigning the final state pions to the intermediate mesons in the 
decay. We find 17% of events have multiple allowable pion combinations. If multiple com- 
binations exist which satisfy the invariant mass constraint on the p° and af (1260) mesons, 
we choose the combination such that the daughter pions of the p° have the larger trans- 
verse momenta. With this selection, the probability of selecting the wrong p and af (1260) 
combination is found to be 27% in the signal MC simulation. 

The dominant background in the B° — >af (1260)7r :F candidates comes from the continuum 
e + e~ — > qq (g = u, d, s, c) processes, since the final states contain only charged pions. To 
suppress the continuum background we apply a Likelihood Ratio (LR) cut. The correlated 
event shape variables are put into a Fisher discriminant 0| to form a single variable, which 
is combined with the cosine of the angle between the B meson flight direction and the beam 
axis, to form a LR. We determine the LR cut value of LR > 0.80 by optimizing the figure-of- 
merit as a function of the LR. From MC studies we find that this cut (LR > 0.80) reduces 
the continuum and b — >c generic decay backgrounds by 97.4% and 76.6%, respectively. 

We estimate the background contribution from rare b—>u decays using large (2500 fb _1 
equivalent) MC samples. Since the B°^af (1320) n T decay also has a four pion final state, a 
significant background contribution is expected from this decay mode. The ai(1260) meson 
has a mass of 1.230 GeV/c 2 and a width of 0.25-0.60 GeV/c 2 , whereas the a2(1320) has a 
mass of 1.318 GeV/c 2 and a width of 0.107 GeV/c 2 |l7J. The upper limit for the branching 
fraction of B° -> af (1320)vr^ is given as 3.0 x 10~ 4 (90% C.L.) 0. From MC studies, 
we find that the 02 (1320) has a large overlap with the ai(1260) in the distribution, 
where 777.3^ is the invariant mass of three pions. However, since the spins of ai(1260) and 
a2(1320) are different, the helicity angle distributions for the ai(1260) and a2(1320) mesons 
differ. A three dimensional (3D) (M hc -AE-cos 9 hc \) fit to the data is used to discriminate 
the signal from the a 2 (1320) background. We note that this method is nearly independent 
of assumptions about the masses and widths of the ai(1260) and a 2 (1320) mesons. 

The signal yields are extracted from a 3D (M bc -A£'-cos^hoi) unbinned maximum likeli- 
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hood (ML) fit to the data. We fit to the candidates in the loose signal region. We define 
the likelihood function as: 

i 

+ f a27T (P a27r (M bci ).P a ^(A^).P a ^(cos^ heH )) 

+ f^u (n^(^b Cl )-n^(A^)-n^(cos0 hcH )) 

+ f(qq+b^c) (^(gg+b^cjl^bcJ-^gg+b-^cjl^-^J-^gg+fe-^^lcOS^hcH))]' (1) 

where i runs over all events in the sample. The coefficients f ain , f a27r , f b ^ u and fr q q +b ^ c \ 
denote the fractions of events from P° — >af (1260)7r T signal, P° — >af(l320)n T background, 
charmless rare b — + u decay background and continuum and generic decay (qq +b — > c) 
background, respectively. 

The P(M hci ), P(AE i ) and P(cos# heli ) are the Probability Density Functions (PDFs) in 
M bc , AE and cos# hel for the signal and background contributions, respectively. The PDF for 
the P° — >af (1260)7r T signal is taken from a MC simulation, with the calibration correction 
obtained from the P° — > D~tt + , D~ — > K + tt~tt~ control sample. The M bc and AE signal 
PDFs are described with a double Gaussian function and the cos# hel signal PDF with a 
third order polynomial function. The PDFs for continuum and b — > c decays are taken 
from the sideband data, where the M hc sideband region is defined as 5.21 GeV/c 2 <M bc < 
5.26 GeV/c 2 , and the AE sideband region is defined as 0.10 GeV < AP < 0.35 GeV. The 
shapes of the M bc , AP and cos# hel distributions for the sideband data are assumed to be the 
same for continuum and generic b—>c backgrounds. This assumption is justified using a MC 
simulation. The M bc , AP and cos# hel PDFs for continuum and generic b —>■ c backgrounds are 
an ARGUS function, a second order polynomial and a third order polynomial, respectively. 
The PDF for P° — ► a^(1320)7r T background is also taken from a MC simulation. The M bc 
and AP shapes are modelled as a single Gaussian and the cos# hel distribution is described 
by a 4th order polynomial function. The b— >u background PDF is taken from MC samples. 

In the fit we float the f ain , f a27T and f b ^ u coefficients and constrain fr q ^ +h ^ c \ according to 

f( q q+b->c) = 1 -(fa 1 7r+fa 2 7r+fb^u)- Figure EJshows the projections of the M bc , AP and cos# hel 
distribution. The M bc plot shows the projection with -0.04 GeV <AP< 0.04 GeV, the AP 
plot shows the projection with 5.27 GeV/c 2 < M bc <5.29 GeV/c 2 and the cos# hel plot shows 
the projection with 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.04 GeV < AP < 0.04 GeV 
(tight signal region). 

A total of 13541 events survive after all the selection requirements and are included in the 
fit. From the fit we find 394 ± 33 P°^a^(1260)7r T signal events, 34 ±38 P°^a^(1320)7r T 
background events and 26 ± 30 rare decay background events. The remaining events in the 
fit are from continuum and generic b— >c background. 

We calculate the branching fraction using the following equation: 



N _ 

£(P -^a±(1260)7r T ) = , (2) 

N BB -e-Br™ b -e PI D 



where N a ± nT . is the number of P° — > af (1260)7r T events obtained from the M bc -AP-cos # hel 

simultaneous fit, N B g is the number of BB pairs, Br sub is the branching fraction of the 
af (1260) decay, e is the reconstruction efficiency estimated from MC, and epm is the cor- 
rection for the difference in particle identification (PID) efficiency for pions between data 
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FIG. 1: M bc , AE and cos 9-fa e i projections for the 3D unbinned ML fit for the selected sample. 
The top, middle and bottom figures show the M hc , AE and cos# hcl distributions, respectively. In 
each plot, the top solid line, the dotted dash line, the dotted line, the bottom solid line and the 
bottom dash line represent the total fit function, the continuum background, the signal yield, the 
B°-^af (1320)7r T background and the rare decay background, respectively. 

and a MC simulation. The estimated values for e and Epiu are 0.061 and 0.962, respectively. 
The branching fraction of af (1260) ^ p°n ± decay is assumed to be 0.5. 

Inserting all these values in Eq. (2) we measure the branching fraction E(B° — > 
af (1260)ttT) = (48.6 ± A.l(stat)) x lO" 6 . 

To measure the ai(1260) properties (mass and width) and the contribution from the 
non-resonant components, we select 1193 B° — > af(1260)7r T candidates with the tight sig- 
nal region to enhance signal events. We fit the ai(1260) mass within the mass window 
0.8 GeV/c 2 <m 3vr < 1.775 GeV/c 2 . For the m 37r fit, we use the following likelihood function: 
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c = II [f^P^M*) + fa 2 ,Pa 2 .(<„) + 



fb->uPb^u( m 3w) + f(qq+b^c)P( 



(qq+b-*c)\ m 3TT i 



(3) 



where i runs over all events. The coefficients f ain , f a27T , fb-* u i anc ^ f(qq + b^c) represent 
fractions of events from signal, B° — > a^(1320)7r T background, fe — > -u background and 
(gg + 6 — > c) background, respectively. The P(m 3n ) are the PDFs for the corresponding 
signal and background contributions. In the present analysis we impose several kinematical 
cuts in the event selection, which produce a non-uniform efficiency in m 37r . We obtain an 
efficiency correction curve from a MC simulation, which is applied to the relativistic Breit- 
Wigner (RBW) function for the ai(1260) and 02 (1320) PDFs. The sideband events and 
MC are used to determine the background shapes for (qq + b — > c) and b —> u, respectively. 
The PDFs for the (qq + b —>■ c) and b — *■ u backgrounds are 5th and 6th order polynomial 
functions, respectively. In the fit we fix the fractions of B° — >af (1260)7r T , B° — >a 2 (1 320) tt t 
and b —>■ u to those obtained from the 3D M hc -AE-cos 8 hel fit |l8j]. We allow the mean and 
width of the a^(1260) RBW function to float in the fit, while those of of (1320) are fixed 
to the world averages • Figure El shows the and p° mass distributions for the events 
that survive after the tight signal region cut. From the fit we determine the mass and width 
of the ai(1260) as 1197 ± M(stat) MeV/c 2 and 305 ± 43(stat) MeV/c 2 , respectively. 




FIG. 2: (a) and (b) p° mass distributions for the events in the tight signal region. In plot 
(a) the top solid line shows the total fit, the dotted line shows the contribution from qq + b — > c, 
the dashed line shows the contribution from B° — > a^ = (1260)-7r =F , the dotted-dashed line shows the 
contribution from B° — > af(1320)7r T and the bottom solid line represents the contribution from 
the b — > u background. 

In order to measure the non-resonant components, we extend the likelihood function of 
Eq. (3) by adding non-resonant contributions from B° — > p°7r + 7r~ and B° — > 7r + 7r~7r + 7r~. 
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The extended likelihood function includes the additional coefficients f p7T7T and f 4w which 
represent the fractions of the non-resonant B° — > p°7r + 7r~ and B° — > 7r + 7r~7r + 7r _ decays, 
respectively. The PDFs for B° — > p°7r + 7r~ and 5 — >7r + 7r~7r + 7r~ decays are taken from MC. 
We perform the m 37r fit in the 0.80 GeV/c 2 < m 3n < 1.775 GeV/c 2 range to extract the 
non-resonant components. In the fit we fix the b — > u and (qq ± 6 — > c) fractions to those 
obtained from the M bc -A.E-cos hel fit. We also fix the values for the mean and width of 
the ai(1260) to 1230 MeV/c 2 and 300 MeV/c 2 , respectively, which lies in the world average 
limits HJ. We apply the constraint: 1 - f b ^ u - f {qq+b ^ c) = / 0lW + f a27T + f pw7T ± / 47r - Thus 
there are four free parameters in the fit: / 0l7r , / a27r , / p7r7r and / 47r . We combine the f p7rn 
and / 47r fractions and estimate the total non-resonant component from B° — ► p 7r + 7r~ and 
5° — ► 7T + 7r _ 7r + 7r _ . We obtain a non-resonant component fraction of 5.3 ± 5.4%, which is 
consistent with 0. Accordingly, we include the non-resonant contribution as a systematic 
error in the branching fraction measurement. 

The measured signal yield from the 3D fit contains several systematic uncertainties. The 
systematic error in the signal yield is estimated by varying each parameter of the fit by 
±lcr from the nominal values. The shifts in the signal yields are then added in quadrature. 
The corresponding uncertainties are listed in Table |I] We assign the uncertainty in the 
track reconstruction efficiency to be ±4.8%. We measure the uncertainty in the efficiency 
of the PID selection of the charged pions to be ±1.4%. A systematic error of —5.3% is 
assigned to a potential contribution from the non-resonant B — > An and B —>■ pmr decays. 
For the selection of the LR cut, a systematic error of ±2.7% and —2.3% is measured. 
We estimate the uncertainty in the number of BB pairs to be ±1.1%. For the efficiency 
calculation with MC we assign a systematic error of ±0.5%. The reconstruction efficiency (e) 
is estimated with MC simulation generated with m ai = 1230 MeV/ c 2 and r ai = 400 MeV/ c 2 . 
We vary the m ai and T ai to our measured values and estimate the efficiency variation, 
and find that the difference is —0.5%. We thus add it as a systematic error. We also 
vary the AE range for the 3D fit to check possible variation of background contribution. 
The AE is varied from -0.10 GeV < AE < 0.10 GeV to -0.30 GeV < AE < 0.30 GeV. 
We find ±4.3% and —1.6% signal yield variation and assign it to the systematic error. 
Finally, we add a ±1.8% systematic error for the fit bias. The total systematic error is 
±7.8% and —8.1%. The branching fraction is measured to be B(B° — > af(1260)7r T ) = 
(48.6 ± A.l{stat) ± 3.9(syst)) x 10~ 6 . 

As a validity check of the procedure in the present analysis we measure the branching 
fraction of the B° — > D~tt + , D~ — > K + 7t~tt~ decay. For this check we perform a M bc -AE 
unbinned maximum likelihood fit to obtain the signal yield. The branching fraction is 
measured to be (2.46 ± 0.05) x 10~ 3 , which is consistent with the world average value of 
(2.76 ±0.25) x 10" 3 0. 

In conclusion, we have performed a new measurement of the branching fraction 
for the decay, B° -> af(1260)7r T , with af (1260) -> p°7r ± and p° -> tt + ij-. Based on 
a 253 fb _1 data sample, the measured branching fraction is B(B° — > af (1260)7r =F ) = 
(48.6 ± A.l(stat) ± 3.9(syst)) x 10 -6 , which is consistent with the BaBar measurement 0. 
Fitting with an efficiency corrected relativistic Breit-Wigner, we measure the parameters of 
the ai(1260) to be m ai = 1197 ± 34 MeV/c 2 and T ai = 305 ± 43 MeV/c 2 . These values 
are consistent with the BaBar measurement and the world average values 
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TABLE I: Systematic errors for the M h -AE-cos hel fit. 



Source of the systematic error 


Relative error (%) 




— a 


Track reconstruction efficiency 


4.8 


-4.8 


PID efficiency 


1.4 


-1.4 


Fraction of the non-resonant components 


0.0 


-5.3 


Continuum suppression cut 


2.7 


-2.3 


Signal shape 


1.3 


-0.9 


B° — ► af (1320)7r =F background shape 


0.4 


-0.3 


Continuum background shape 


1.7 


-1.5 


Rare decay background shape 


0.4 


-0.4 


Number of BB 


1.1 


-1.1 


MC statistics 


0.5 


-0.5 


m ai parameters in MC 


0.0 


-0.5 


AE cut selection 


4.3 


-1.6 


Fit bias 


1.8 


-0.0 


Total 


7.8 


-8.1 
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